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CHAPTER-I 


1, INTRODUCTION 

The study of 'Fluid Dynamics' is directed to the 
behaviour of a fluid in motion. The liquid and gas states 
are referred to generally as fluids. Some of the 
examples of application are : flow of water through pipes? 
motion of an aircrafts or a missiles in the atmosphere. The 
study also yields the methods and devices for the measure- 
ment of various parameters, e.g. the pressure and velocity 
in a fluid at rest or in motion. 

The first notable works on fluids appeared in the 
seventeenth century. In 1687, Newton in his book 'Principia; 
deals with the influx of a fluid, the resistance of a fluid 
and the resistance of projected bodies and the wave motion. 
Daniel Bernoulli, in 1738 determined the relation between 
pressure and velocity and formalized this in his theorem. 

In 1743 his son John Bernoulli applied the momentum principle 
to infinitesimal elements. Euler realized the significance 
©f the Bernoullis' work and used it as the basis, setting 
down the fundamental equations of motion and equation of 
continuity for an ideal inviscid fluid in 1755. An alter- 
nate form of these equations was given by Lagrange# in 
1781 and 1789. Equatiohs including viscosity were derived 
by Navier in 1822 and Stokes in 1845. 

In 1858 Helmholtz published a paper on vortex motion 
and in 1868 another ©n free stream line potential flow. 

Near the end of the nineteenth century, fluid flow 


: ; ^ 

itself began to be extensively observed and investigated. 

The modern trend in fluid dynamics is to investigate the 
flow of electrically conducting fluid at very high temper- 
ature. 

At a microscopic level, the three states of a 
given s\ib stance are different because of the difference 
in the intermolecular distance. In many cases at normal 
conditions the molecular distance of a fluid are less than 
the minutest of any physical dimensions of practical interest 
As a result, we interested in the statistical average prop- 
erties and the behaviour of large numbers of molecules, and 
not in that of individual molecules. That is, macroscopic 
and not microscopic, properties are ®f interest. In Fluid 
dynamics as individual molecules are not being considered, 
the fluid can be regarded as a continuous medium and so the 
physical quantities such as mass, momentum etc. of the fluid 
contained in a very small volume are regarded as being 
spread uniformly throughout that volume. 

In dealing with the gases at very low pressure, as 
in the upper atmosphere, or at very high temperatures such 
as in a plasma, the continuum concept of fluid dynamics 
must be violated and the study that has to be based on the 
behaviour of individual molecules (i.e. on the macroscopic 
approach) [ i] . ^ ^ ^ ^ 

In our studies we frequently refer to a ' small 
element of fluid* which is always supposed so large that 
it still contains a very large number of molecules, as 
fluid is a continuous medium. So when we take of an i 
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infinitesimal element of volume we always mean that which 
are physically infinite small. Such an element is called 
a fluid particle. 

In continuum dynamics we as sum that the macroscopic 
fluid properties, for example mean density, mean pressure, 
vary continuously with (a) the size of element of fluid 
considered (b) the position in the fluid, and (c) the time. 
In (a), the variation becomes negligible as the element is 
physically very small. Thus, fluid properties density, 
pressure and velocity are expressed as continuous functions 
of position and time only. On this basis, it is possible 
to establish equations governing the motion of a fluid, 
which are independent of the nature of the particle struc- 
ture. So gases and liquids may be treated together. 

There are two distinct methods of specifying the 
flow field£ 2] "* 

(i) LAG RANG IAN METHOD ; In this method, the flow varia- 

bles (velocity, pressure and density) of a selected fluid 
element or particle are described. If is the position 
of the center of mass of the flviid element at time t^ then 
the basic flow quantity in the Lagrangian description is 
the velocity ▼ ( r^, t ) . This method is also referred 
as 'individual time rate of change' . 

(li) EULERIAN METHOD ; In this method, the flow quanti- 
ties are described at all points of- space occupied by the 
fluid at all times i.e. flow quantities are defined as the 
function of position in space (r) and time (t) . The basic 
flow quantity is the vector velocity t (r, t). This method 
corresponds to ' local time rate of change *. 
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with some property of the fluid (it could be density# velo- 
city etc, )• Suppose fluid particle has the position 
PC x#y#z) or P(r)at time t. Keeping this point fixed the 
change in f is during the interval of time ^;t is, 
f( x,y,z,t + j t) - f( x,y,z,t). 

Hence theMocal time rate of change' is given by, 

. 


^t 


lim f(x, V, z,t-t- St) - f(x, y, z,t) 

i5t-^o ' St 

As the point P is fixed the local time differential 
operator is not carried along by the moving fluid. 

Now let at time t+st the fluid particle, which was 
at the position (x,y,z) originally, is in the position 
(x + ust, y + vjt , 2 + wst) where u,v, and w be the velocity 
components at the position at time t. The corresponding 
change of f is given by, 

f(x + uit, y + ust, z + Wit, t +5t) - f(x, y, z,t) 

©r f (r + sr, t >fc s t) - f (r, t) 
and rate of change is. 


df 

It 


lim ■ 


fCr + &r, t 4-it) - f(r, t) 


This gives the ' individual time rate of change*. As the 
point P is moving it gives the rate of change which is 
carried along by the moving fluid. 

Since, 


f “ f(x,y, z,t) 

d£ . ^ ax + dy + If 


dz + — dt 
dt 


% / . 


SO, 


df _ 3f af . Bf 

dt *" at ax dt ay 


+ if- 

dt 32 


dz 

dt 


If V = [u, v,wj be the velocity of the fluid 


dx 

particle at P and = u. 


dt 


dt 


= V, 



first term qn the right hand side represents the lo< 
rhanae of t. and the second term the convective 


If V is introduced for f in the above equation -cm 
al' derivative of velocity with respect to time is 


It is this term xn the Euxer rex< 
which is responsible for the non 
ions of motion of fluid dynamics 



2. FUNDAMENTAL EQUATIONS 

In order to study the details of the fluid flow, we 
wish to find the density distribution, pressure distribution, 
velocity distribution, states of the fluid, etc., at all points 
of space occupied by the fluid at all times. Hence a knowle- 
dge of three velocity components (u,v,w) , the temperature T, 
the pressure p, the density -f , etc., of the fluid which are 
functions of position in space (r) i.e. (x,y,z) and the time 
t, is needed. We obtain relations connecting these unknowns 
and which would explain the particular problem of the fluid 
motion. We call these relations as fundamental equations. 

The flow of a compressible, non-heat conducting and 
inviscid perfect fluid is described by the four physical 
variables pressure p, density f , vector velocity v, and 
teraperature T. So there must be four fundamental equations 
to find these variables. These are^3 ], 

Equation of continuity, 

+ div (-f v ) =0 


it* v; 

Equations of motion. 


0 in one dimension 


M - it* ‘ = - T 

When the external force is negligible, 
Energy equation. 


PS- . Ph , Js. J2JE 

Dt Dt f Dt 
for the adiabatic flow. 





and equation of state for perfect gas, 

P = f RT 

when there is no intermolecular attractions. 

NOW we consider the behaviour of the fluid flow at 
high temperature, «here some of the fluids become conducting 
and the interaction between motions of the conducting fluids 
and variations in electromagnetic fields may not be neglig- 
ible. To describe this interaction phenomena of the flow 
field with the electromagnetic field we .find the velocity 

vector ▼# pressure p, density f , temperature T, and the 
magnetic field vector H, which are the functions of the 
spatial coordinates (x,y,z) and the time t. To calculate 
these unknowns we must find the relations which are the 
fundamental equations for a conducting compressible fluid. 

In these equations there are coupling terms between the 
electromagnetic aiid fluid dynamics phenomena. We shall 
derive and discuss these equations! 4 J , 5j > [® ] ^ 

equation of continuity t 

This equation simply expresses the law of conserva- 
tion of mass. • 

The quantity of fluid entering a certain volume in 

space must be balanced by that quantity leaving i. e. matter 
is neither create# nor destroyecie 

We now formulate this principle mathematically. 


i: 


Let V be any arbitrary volume fixed in space, f 
bounded by a surface S, and containing a fluid of density 

f o The volume element (JV is small so that f can be regar- 










% 
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ded as constant through it 


The mass of the fluid with in the volume V is j j ^ 
The rate of generation of the fluid with in the volume is 


is constant with respect to time. 

If the volume V occupied by a moving fluid, the fluid 
enters V through parts of its boundary surface S and leaves 
through another part* 

Let n be a unit outward normal vector drawn on the 
c>Tn<=.«i- as. T*Vi« normal velocity is The total 


Using Gauss' s theorem. 


The sum of the net outward convection of mass plus 
:e of Generation, of the fluid with in the volume must 




Since this is tme for arbitrary elementary volumes. 


( 1 . 1 ) 


The Maxwell* s electromagnetic equations for a conduc 


ting medium are. 


curl E = - ^ where B = H 

— - 

curl H « 4 ‘H’ J + 

div B =0 

and, D » 4 IT q 

Differentiating the fourth Maxwell's equation 

div C^) “ 

using the second Maxwell* s equation it becomes, 
* div ( curl 1 - 4 -n j ) = 4fT 

Third Maxwell* s equation is 


we have. 


or + div 3 = 0 

or, 

If the charged particle moves with velocity v 




EQUATION OF STATES 



The electrodynamical state relation is a simple relation 
between the current density J, fields and fluid motion. Or, 
the Ohm' s law is, 

J = -- [e + Av X 1 ] + q V 

Where <r- is the electrical conductivity and qv is the , 
current depending on the motion of the net charge q. 

The state of a compressible fluid is defined by the 
pressure p, the entsopy S, the internal energy £ , the absolute 
temperature T, the mass density -f and the specific volume Vv 

An equation of state is formed by expressing any one 
of the thermodynamic variables in terms of the other two quan 
tities mentioned above. 

For example, 

p = p(f , T) 
or, E - E(P, -f ) 

Pqj- a perfect fluid Cno intermolecular intractions) , 

having negligible viscosity the equation is, 

p = R f T 

Where, R is the universal gas constant. 

Assuming the process to be adiabatic and isentropic we have 

■■ £ ■= c T ' 

V 

£ « ^ as p = -f RT 

Also, Cp - = R 

or. 



This is the caloric equation of state of the medium. is 
internal energy per unit mass* 

When a fluid element changes its state isentropically 
the state equation is. 


EQUATIONsOF MOTION! 


This equation express the laws of conservation of 
momentum in the fluid, Te obtain equations deriving this we 
proceed as followings - 

We consider a fluid mass in motion at time t occup- 
ying a volume V and bounded by a surface S. Let S V be a 
small element of volxame. If the fluid is of density -f , 
mass of the element is f 5 V and it moves with velocity ▼(r. 
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The inertial force on the element is f VV( • This 

force equals the rate of change of linear momentum of the 
element. For the whole fluid mass. 


m 


f iV 


force. 


The force on the body is the sum of body and pressure 


Consider a surface element 5S and that n be the unit 


outward normal vector to JS, Let p is the fluid pressure at 
any point on the element ^S. Then the force on it due to the 


fluid outside the particle is, 

p( - n) SS 


(For pressure acts along inward) 


Hence the total pressure force is, 

- J Jnp dS 

using Gauss's theorem. 


P 


- Ill 


grad p dV 


When the fluid mpves in an electric and magnetic field, 
the body force F„ consists of three parts: gravitational, electric 
and magnetic. 

The gravitational body force on the element is f fiV 
where g is the gravitational acceleration. 

The electric body force on it due to an electric field 5 
of intensity E would be E q SV, 

The current flowing through the element is, i = J 5 A, - 
Where S A be normal cross - section of fluid element. 

Hence using Biot - Savart law, magnetic force in the 


element is. 




i il X B 
? 2A Si X S 








r 






» ( J X B ) 

» ( J X H ) SV 

Hence the total body force Fg on the fluid mass is, 
Fg « I j J ( -f f + 7^ J X H + qS ) dV 


Now, 


** F_ + F 

I B . p 


j I" . f^ + ( f g + XH + ql ) ~ gradp j dV = 0 

' ' v 

The above equation is true for arbitrary elementary volumes 


and so. 


f ^ * - grad p + -f f + /*“ J X H + qJS 


This is the general equation of motion for a conducting 
fluid. Where Euleri an derivative is, 

SB -iU_ + ( ▼• ^ ) 

The term involving E can be neglected, as compare to 
the term involving magnetic field H • As the order of magni- 
tude of qE is very-very small as compare to the order of 
magnitude of ^ J X H i.e* 

0 ( 1 qE j) « 0(1 M J X 1 1) 

So the modified form of equation of motion is, ■ ■ ; 


We know. 


» - grad p + -f g + ? X H 

curl H “ 4 fT J + (: 














The term fe can be neglected as compare to curl H 

dt. 


0 ( 1 t ^ 1 ) « 0 ( 1 Curl H 1 


Curl H * 4Tr J 

Then the magnetic force is, 

^ j X H = Curl i X H 


H X Curl H 


Since, C using vector identity ) 


H X Curl H * grad (— y-) 


(H. ) ii 


Then, 


>^J X H * - grad ( 


^ ^ Cii.^ ) H 


So equations of motion can be written as. 


^ i? = - grad p + -f f - grad ( 

In summation convention form, it can be written as. 


''j > + ''j -fx^< ^ 




where 


i is the dummy index 
j is the free -index 
i, j = 1, 2, 3 
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EQUATION FOR THE VARIATION OF THE MAGNETIC FIELD 


This equation expresses how the magnetic field vary or 
the time dependence of the magnetic field. 

The Maxwell's electromagnetic field equations for a 
conducting mediiim are, 

div E = 4Tr q 
div H = 0 


Curl E = - 


1 3H 
o at 


and. 


.. 47r — 

Curl H » ■— J + 
G 


The first two equations are derived by applying Gauss's 
theorem to a closed surface. Third equation is the mathematical 
formulation of the Faraday's law of induction and last equation 
describes how the magnetic field ii depends on the conduction 
current J and the displacement current • 


Since, 


curia = if 3 + 


1 S E 

By neglecting the term ~ above equation reduces to. 


•w* 'IT mm 

Curl H » 

e 


Taking Curl of this. 


Curl Curl H *= Curl J 

c 


grad div H - H « -^~Curl fo-(E + 


where. 

Since, 


J*= cr(E + -ivXH) 
c 


div H ^ 0 










- V H » 


\lT<r- r 

C L 




Curl E + -i Curl (v X 5 )J 


; ". 5 ' ^ 

f I*' ■*' , V 




rV- '''fir'-;:." ' 


Also, 


2 — 

S7 H 


Curl E *= - 

4 TTo- f 

“5T- I - 


at 

we may write. 


« Curl X H ) + 


c 

at 

SH 


at 

+ 

>.2 

c 


4170" 


2 



M 


+ Curl ( T X 5 ) 


where ^ 

• AlT<r 


« (H. V )▼ - (v. v ) i - H ( V*V ) + 1 H 


In summation convention form, it can be written as, 

■it ‘«J> = <«i > ’'j - <’^1 ^ ) Hj - v^) + 1 ( 4?>“l 

^ ■ . i X X 

■ ■ ■ 

This is the magnetic field equation for a conducting compre- 
ssible fluid. 

For' a fluid at rest it reduces to the diffusion equation, 

-|S- » v| ^ 2 if 

*1 is called magnetic diffusivity. For an equation of 

this type the rate of decay of the magnetic field is very rapid. 

And when the conductivity is so large that the term 
2 — 

^ H can be neglected, the behaviour of the magnetic field 
is given by. 









IIIM 



= Cmrl ( ▼ X 1 ) 


» ( H. ) V - ( T. ) H - H ( V 


In this case the rate of decay of magnetic field is much 
less therefore this is the case of high conductivity. 


EQUATION OF ENERGY: 


The law of conservation of energy leads to another 
equation of fluid flow - the energy equation. 

The law of conservation of energy is equivalent to the 
first law of thermodynamics, which states that, if a small 
quantity of heat is added to a simple system it is used up in 
changing the internal energy of the system and in the external 
work done by the system. Mathematically this law can be expre- 


ssed as. 


d£ o dQ . pci( ) 


Where dQ is the heat added per unit mass, d£ is the 

increase in internal energy per unit mass and pd(--~) is the 

1 

amount of work done during the change of volume d(— • The 
heat dQ is obtained by heat conduction. 

For a perfect gas, the molecules have very negligible 
volume and there are no mutual attractions between the indivi- 
dual molecules and hence no potential energy. 

Then energy per unit volxame of fluid flow is the sum of 
the kinetic energy density due to translatory motion of the 

. ' ■ ’I ■ 2 ■ ' 

molecules, —2“ f v and the internal energy density f£ . 


I 




■ ■■ , 




The equation of motion is reduced to. 


f = - grad p + qE + J X H t neglecting the gravi- 

tational force. 

Taking the scalar product of the above equation and the 
velocity vector V we get for the kinetic energy, 

_v 2 

f « , (y, V ) p+q (E.v"") +/A[r. (JXH)] 

The internal energy of the perfect gas depends on its 

f 

temperature T only, 
we know. 


dt *= 

dQ - 

dt 

dQ 

lt~ “ 

dt 

d£ 

dQ 

■III]*— Biwii 3CS 

dt 

dt 


c-f, 


f 2 dt 


The equation of Continuity for a compressible fluid is^ 


+ f div V = 0 


d£ _ dQ p 


V . V 


It is an experimental fact that heat flow across an 
element of fluid surface in Jt time is proportional to t: 
gradient in the temperature i.e. 

SQ oC - is St -l-i 


■ ..ft*' ■■■I 


t';-' ,* 

r.' I ' - . ' . . 





7 '"j 


' , ; ■■ -I ■ r' ■■ ‘ '■ 


liiur >.'i 




§Q » - k 5A St 

• 

where k is the coefficient of thermal conductivity. The 
negative sign implies that heat flows from points of higher 
temperatures to points of lower temperatures. 

Consider a closed surface S enclosing an arbitrary volume 
V, The heat flow out in %t time is, 

” Ij ^ * - St Jj k( V T). dS 

= - St liJ K7>(k VT) ciV 
The heat added to the system per unit volume is then. 


f SQ a St V . [k'7T] 


Therefore, as Jt 


So, for the internal energy 


t - P C ) 


f [ k v tJ - p ( V . v ) 

Hence the sum of kinetic and internal energy gives. 


( £ 4 ^ -(▼.v)p-.p(v,T)+V', [kv t] + 

+ q(E. V ) [ ▼. ( J X H )] 


( £ + ^) 


5 -) « -V.(pv) + V,[kVTj + q(I.v) +>^[v. ( JxR )J 













Now we find the ordinary energy equation of an electromagnetic 


field. 


When the fluid moves in an electric and magnetic field 


it experiences a force known as the Lorentz force and experi- 
mentally its value is, 

P = qi + •— ( J X B ) 

■ c ■ 

Where qE is the electric force per unit volxme due 

to an electric field and ( J X B ) is the magnetic force per 

■ ' / 

unit volvirae due to magnetic field. 

The rate of work done by the Lorentz force P is. 


P.v » q ( E.v ) + (J X B)J 


(1.3) 


This force is the cause of the motion of charges if 
they are free in the field of E and B* If r is the velocity 
of motion of the average charge density q at any point, we have 


« i. J + ^ [( V. ( ▼ XH )] 


p.v » E. J as V. (▼ X 5) = 0 


(1.4) 


Now we transform the right hand side of above equation 
(1,4) into a function of E, 5, B and H with the help of Maxwell 
equations. ^ 


Since, 


Curl H *» ^-2— J + -i — ~ 

CJ c d t 


:. Curl H « 


(I. 5) t ~ t E. 

c 










> W,, , 


Using the vector identity, we have 



E. Curl H « H. Curl E - div ( E X 5 ) 
with the aid of this result 


E.J «= [h. Curl E - div ( 1 X H )J - (1. 

, Since, Curl E ■ - ~ 

c at 

E.J = - div ( E X H) - [ CH. -it ) + ( E. -|| )] 

If the medium is isotropic then 


We write. 


Since a part of the power density is the rate of vari 
2 2 / 

tion of ( 6E + /*-H ) Y , this quantity must be the sim 

of energy densities in electric and magnetic fields. 


Using (1.3) and (1,4) equation (1,2) becomes 





1 

■ 


Adding equa-tions ( 1.5) and (1.6) we get 


f -iL. / £ + 
1 Dt ' ^ ^ 


-»• 


* - V, (v p) + V. (k V T) - V. (E X H) 

This is an equation for total energy per unit voliame of 
the flow field, when there is no heating by radiant energy. 

The terms on the left hand side of above equation give 
the rate of change of energy per unit volume in the field. The 
right hand side terms of the above equation have the form of 
divergence so that, on integrating through a certain volume, 
all of them reduce to surface integrals, . 

The first term represents the work done by the pressure 
force, second term represents the heat loss through conduction 
from the surface and the last term represents the out flow of 
electromagnetic energy through the surface where E X H is the 
Poynting vector. 

The above equation of energy can be written as. 


3 , . 1 ^ .2 . . 3 / _iE“ . M X . 

3 1 ^ ' V 2 ^ ^ ) + ^t ( Qjf * 8 TT ^ 

+ div |^ fv ( t + "I- ) + prj = div(k grad t)- -^divCE X H) 

For adiabatic flow, no heat is added, conducted or 
radiated from the flow field, i.e, 

dQ * 0 

The energy equation for the adiabatic flow field. 
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-gHOCK WA VES, their EXISTgung 

If a small disturbance ( l.e. v,lth an infinitesimal 
amplitude and small velocity) is created within a non-vis 
laentropic compressible fluid it will propagate throughout 
fluid as a wave motion and with the velocity of sound 
to the fluid Without suffering any distortion. 

in this case we get the wave equations, which shows 
both the density and velocity variations fellow 
patterns, derived from the linearizing of th 
fluid motion. The solution f (x. t) and u 
single valued 


In sxiramatxon convention form. 


) + div C i v( £, + 


fe + 


Where, 


snd. 


Mjk 







in Euler* s equation of motiqn may be reglected. 

If the assumptions of infinitesimal amplitudes 
and gradients are removed# then the wave velocity will not 
be a constant and a simple wave velocity will distort as 
it propagates. ' 

In this case as velocity is not small# we have to 
solve the complete non-linear equations of fluid motion. 

The method of characteristics can be applied to solve such 
hyperbolic type equations. 

A simple analytic solution of these equations of 
fluid motion may be find such that the density is a function 
of velocity only. After some manipulation we get equations 
in u and f , General solutions of these equations are# 


u - f^ [ X - (u + a) tj 

f - £2 [ X - (u + a) tJ 
where f^ and f 2 are arbitrary functions. 

We consider the solution# 

u = f^ j^x-(u+a) tj 

This equation shows that the disturbance is propagated 
at an instantaneous velocity u + a# instantaneous because 
this velocity is a function of time. If the increment in 
velocity u is much smaller then speed of sound# then the 
solution is the case of sound wave i.e, the linear case 
and the curve u = f^ ( x - at ) does not change its shape 
as the disturbance propagates. But if u is not small the 
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shape of the wave is distorted as it propagates. 

The resulting distortions in velocity distribution 
are shown in figure. 


FORMATION OF SHOCK 


At a time t^< the points of high velocity move faster 
to the right than those of lower velocity. Thus the creast 
has moved faster to the right than the trough and so the 
profile is disturbed. The compressive part of the wave, 
where the propagation velocity is a decreasing function of 
X, distorts to give a triple valued solution for u(x, t), 
which is physically impossible for longitudinal wave. 

Actually, the difficulty is over come by the formation 
of a * shock wave * • In such cases, a large variation in 
pressure and density occurs in a very narrow region in which 
the flow variables change rapidly and the fluid no longer 
undergoes isentropic changes. The thickness of this shock 


wave region is very small and so we may considered the 
shock wave as a surface of discontinuity for many practical 
problems of inviscld fluids/ [ 4] » [?] ). 

The propagation of shock is faster than sound when - 


observed from one side of the discontinuity and less than 
that of sound when observed from the other side. Hence 


the velocity of the shock will be supersonic viewed from 
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ahead and it will be subsonic viewed from behind. 

The formation of these shock waves causes a great j 
noise equal to that of an explosion# usually termed as 
supersonic bang so the effect of shock must be taken into 
account in the design of aeroplanes, pipe flow, super- 
sonic flight of projectiles and so on. 

For the flow of an inviscid and non-conducting gas 
the laws of conservation of mass,mementtan and energy are 
originally formulated in the differential equation form 
as it IS assumed that the flow variables defining the 
flow are continuous functions. Flows are also possible, 
however, for which discontinuities in the distribution 
of these flow variables accur. So the conservation laws 
can also be applied to such discontinuous flow and hence 
across a shoeko 

The condition for the existence of shock waves 
which may be called the jump conditions, relate the velo- 
city, pressure, density and temperature in front of the 
shock to those behind of them. The jiimp conditions are 
the simple consequence of the laws of conservation of mass, 
momentum and energy across the surface of discontinuity 
and the equation of state of the medium through which the 
shock is moving. 


4. JUMP CONDITIONS ACROSS A DISCONTINUITY 




A differential equation can be represented 


(1.7) 


r 








4 * 



per unit time 


Of some state of medium 
and K are continuous and differentiable. 

Let x,{t) and xj(t) are continuous and diffe: 
ion of t and xft) < x^(t) for every t. 

Integrating (1,7), 

oCiCt) 

J at - F(x, ) = - I K 

Now we know that Newton-Leibnitz formula is, 

t) dx = [ dx 4- H (x^/t) 4^ - H( 


Using this formula, (1#8) becomes 


Where P(x^) is the outgoing flux and 

incoming "flux.'. 


amount of state of the medium 

L.H,S. H dx) is equal 

'*•1 

X and incoming flux, and 


some contribution of 




H « H(U) , F » P(U) 

then (1.9) represents the conservation of mass. 


H dx = F{x.) - p(x, ) 


and equation (1.7) is in the conservation fornix 


The most general conservation form for differential 
equation is. 




vi-t I ^ j ~ 0 , n 


1 , 2 ,- - - ( 1 . 10 ) 


Suppose there is a discontinuity at x = X(t) and x, 
and Xa_ are chosen so that x, < x(t) < x^^ . Suppose and 

VI ■ 

U ana their first derivatives are continuous in X(t) > x 

and in ^ x > X(t) and have finite limits as x — > X(t) from 
below and above. 

Integrating (1.10), 


■k ^ 


■iTT [u (x,) - u (Xj,) j 


Also using the property of line integral, 

xltj ^ 

4^ dx = f dx + f 


( 1 . 11 ) 


Now we use the Newton-Leibnitz formula on the left hand side 

7^ . "w o.. 


ol'<. 


T^.4, 

u-’d-. + jj uV -u'fx-:;*;; ^ 


N. ' . 














Where U (X , t) , u" (x"*” , t) are the value of U** (x,t) as 
X — > X(t) from left and right, respectively, 

^4^ = u" dx t u’^dx + S ( ) (1.12) 

"*■1 '*-1 xLt) 

where, ^ ^ 


U (X “ , t) 


'll '^y ’4" % ^ 

U, , U (X ,t) — ^ Ua 


On comparing (1.11) and (1.12) we get, 
xtt) , V, r A 

u-'d- ^ >d-d- - s (uj- = Jrr " 

at ^ J ‘ 


Since is bounded in each of the intervals separately 

the integral tends to zero in the limit as x, — > X and 


x^— > X 


Therefore, 


s(u;-ujj = iTPi ( 

= ^ ( u-'tx-J - olVj] 

, *1 H . Vi / '^+’f H4- f % 

This condition may also be written as, 

- S[U^J ^ =» 


where the brackets [ ] indicate the jump in quantity and . 
^ is the speed of propagation of the discontinuity or shock 
front. 
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This gives the location of the discontinuity. So if 


the cunservaticn equation is, 

5H . ^ _aL =: o 


at a'>*- 


(1.13) 


then the jump relation is,[ 

- 5 1«] [f] 


a. JUMP CONDITIONS IN ORDINARY FLUIDS 


(1.14) 


Consider a stationary shock front (discontinuity) 
separating two regions* i.e. ^ — ^ so, 

the jump condition (1.14) becomes. 


(1.15) 


■ The equation representing the conservation of mass, 
momentum and energy across a surface of discontinuity are 
given below. 

Equation of continuity is, 

l|_ _i— (^-f i>;J 

at V 

which is in the form (1.13) . So the corresponding jump 


condition (1,15) is. 


[ f = o 




i.e. ~ * 


(1.16) a 


Equation of motion is 












f ^ j + i ~ 


^ t-fiJ-jj ‘^3 ■*" ^ 


(using equation of continuity) 


^ (-fi>jj -*- -1$^ -^- f 5^- ^3 =• - 

or 

4 WJ 

which is in the conservation form, jiomp condition is 

[ ( -f l>i l> j t- l»j VI, J o 

using ( 1 . 16 ) a we have. 


[^]^i 


Equation of energy is. 


( 1 . 16 )b 


I <ll» 

f dt 


In conservation form it can be written as 



Where h + 


£ = ^o- ^ . 
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Jump condition is, 
f t 

I I ftr-ij 






If ■?) '^i ~ then the constant factor 

may be dropped in the above condition and we have. 


f - 

L fu-) 


r] - 


(1.16) c 


The equations (1.16) are the jump conditions (Rankine- 
Hugoniot) when the shock is at rest[8j,The subscript 1 and 2 
denote a variable ahead and behind of the shock front, 
are the components of the unit normal to the shock and bracket 
denotes the difference of values in the two sides of the shock 
surface of the quantity enclosed. 

b. JUMP COKDITIQNS IN MAGNETOGASDYNAMICS 

In an electrically conducting fluids in the presence of 
magnetic fields discontinuity (i.e. a shock wave) in flow vari- 
ables can exists. The study of magnetohydrodynamic shock 
waves was begun in 1950 with the paper of P. de Hoffmann and 
Teller o] . Since then continued interest inspired by 
astrophysics, by flight at the outer edges of the atmosphere. 




etc. has produced many papers describing shock wave properties. 
The basic properties of magnetogasdynamic shock waves are deter- 
mined by the conservation laws. In the presence of magnetic 
field the relations connecting the flow and field quantities 
on the two sides of the shock surface are as follow [loj , 


Equation of continuity is. 


which is in the fom (1.13) so the corresponding jump condl- 


tion (1.15) is. 


[ f - 0 


[fa^J=o 
i.e. =. 


(1. 17) a 


The equation of motion is. 


since 


-t InT C") AiJ 




i“ Hr 


. 3 u. ==. - Hi ^ 

UTiCl ^l A-AiA ‘ V J ^'-■1 




as H = D 


Using these results equation of motion can be written as. 
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which is in the conservation form. So the jump condition is. 




-f 


U* t^vf 


-4-- HH, 


,J=o 


l^-f U- livlj +■ ~ 




(1.17)b 


Also we know 


1. e. 


c^n^ H =i » 


_ 2 _ H.' = 
i-t-i ‘ 


Jump condition is. 


Lh.-'’*J 


[ H.J = >= 


(1.17) c 


Equation for the variation of the magnetic field (when <r--joo) 
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jump condition is. 


. - HiU-jj =: O 


- «i =» 

[ ff D-J = o 


(1.17)d 


The ecjuation of energy is^ 

^ [fUi- I |>(P+ n+ 

Ohm* s law for a moving medium is, 
J = -h ^ dxifj 

or _ .. _ ■_ . 

E t=. J — l?- X H 


■^ixH = ^ H x(i>XHj - ^ H 


also we know 


CuJtl H = J- 


(neglecting the displacement 
current) 


^ EkH = HXtUXHj - Mx(vxHj 

in the case of high conductivity it reduces to 

^ XUXHj 

or ■ 

using this result, energy equation can be written as, 

^ I f -^jj t- ^ (>(^*7- -ffc - (iTiq; ijjj =. 


** / 


























jTjmp condition is. 


1>* c v* + -f £ ■»■ 2- -f - iSW ^ "'i ■ J j “i'J 




(1.17) e 


The equations (1.17) are the shock conditions in the 
frame of reference in which shock is at rest. The equations 
(1.17) are also called as the MHD Rankine-Hugoniot relation, 
which can be written as, from equations (1.17) 

K ~ sttK [(T'+'jV (’"•'J’ftj 


ll+t) lo, i- (r-tjhj. 

t- CY'*-') 


aiT-fi [(y-'jh t(Y+i)i=»4j 


For H =o this relation gives the result for a non- 
conducting fluid. 

The usual situation is that the flow ahead of the shock 
is known and these conditions are used to determine the flow 
behind or to determine the flow quantities in terms of one 
of the flow quantity behind. 













5. radiation PHENOMRKrnw 


At very high temperature, radiation can be considered 
as a continuous emission of energy in the form of electro- 
magnetic wave which propagate in vacuum with the velocity ' 
of light. This energy is called as radiant energy or ther- 
mal radiation. Whereas, according to • quantum theory ‘ 
the radiant energy emitted or absorbed is not continuous 
permitting all possible values, as demanded by the wave 

theory, but in a discrete quantified form, as integral 
multiples of an elementary quantum of energy, photon or 
light quanta. The amount of energy in each quantum being 
given by the product h y , where h is Planck's constant 

and V the frequency of the radiation. 

Thus the quantum theory proposes the particle charac- 
teristics of radiation, while the classical theory the wave 
characteristics, both being required to understand the 
complex behaviour of radiation. It has also been recognised 
that radiation reveals itself in different types such as 
the electrical (radio) waves, infrared, visible, ultravio- 
let, X-rays and gamma rays. This theory of thermal radia- 
tion can be applied to understand the processes which take' 
place in stellar media, to explain the observed luminosity 

of stars and nuclear explosions and also to high temperature 
fluid flow. 

Radiative transfer and radiative heat exchange have 
an influence on both the state and the motion of the fluid. 
This influence is caused by the fact that fluid loses or , 
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gains energy by emitting or absorbing heat. The state of 
the fluid can be described by the fundamental equations 
which, in the presence of radiation field, must include 
the interaction between the radiation and the fluid. There 
are three different thermal radiation effects on the flow 

fieid([i:0> [lilj iial ) 


lergy Dens 


Er » 


The radiant energy density per unit mass of the fluid 


is given by. 


■! 


I 

I ' 
■ 


Where a is Stefan-Boltzraan constant. 


Radiant Pressure Pj.* 

According to the well known result of classical 
electrodynamics the pressure of a radiation field is equal 
to the one- third of the radiant energy density i.e. 


iiM 


I 

1 

li® 

<1 f'' 


' 4 - 

P » ny E_ 

- .i 

The radiant energy density and radiant pressure 
become comparable with the energy density and pressure of 
the fluid only at extremely high temperature or extremely 
low gas densities, , 

Radiation Flux F_ i 

' i ^ - 

The net amount of radiant energy passing through the 
surface per unit area per unit time is called the radiant 


mim 











in the energy equation the radiant energy density 
must be added to the energy density of the fluid and 
also this equation requires the introduction of a term 
describing the energy losses by radiation.The energy equa' 
tion (for adiabatic flow) then becomes. 


6. SPHERICAL AND CYLINDRICAL SHOCK WAVE 


Consider the propagation of a shock wave, through a 
perfect gas, of great intensity (i.e, very strong) resulting 
from a strong explosion i.e. from the instantaneous release 
of a large quantity of energy. 

When the energy is suddenly released, in an infinitely 
concentrated form and distribution of density, pressure etc. 
depends only on the distance from some point then this is 
the case of spherical shock wave. 


When the energy is suddenly released along a line and 
distribution of all quantities is homogeneous in some direc' 
tion and has complete axial symmetry about that direction, 
then this is the case of cylindrical shock. /. 

Since we consider the symmetrical flow (centrally or 







r and t* So the one*»dimensional fundamental equations 
governing the adiabatic flow areffsl, [ 14 ]) 


Where j = 0,1# or 2 stands for plane, cylindrical or 
spherical case, respectively, and r is the radial distance 
from the centre in spherical case, the radial distance from 
the line of explosion, in the cylindrical case and x coordi- 
nate in plane case; u is radial velocity; f , the density. 

Similarity principle may be used to reduce these 
equations to ordinary differential equations. 


7. SIMILARITY PRINCIPLE AND SELF-SIMILAR GAS MOTION 


It is not always possible to solve non-linear differ- 
ential equations describing the physics of a motion or a 
process by using mathematical techniques. An approximate 
solution of such a problem can be obtained by solving a 

I- , -'i 

similar problem which is easier to solve. 

We consider the problem of one-dimensional adiabatic 
flows of a perfect gas with constant specific heat, with 



(, 1 . 10 ) to (1. 20). These gas dynamic equations contain 
five dimensional quantities p, f . u, r and t. The • 
dimensions of three of which are independent; for example 
density# distance and time, so the equations admit three 
independent similarity transformation groups of quantities 
Dimensional analysis can be used to obtain these groupings 
By the successive application of thebe three groups we can 
obtain solutions for the different flow similar to each 


tial and boundary conditions of the problem. 

The motion Itself may be described by the most general 


and u(r,t)o These functions also contain the parameters 
entering the initial and boundary conditions of the problem 
They do not depend upon the position r and time t independ- 
ently but are functions only of the combination ( r/t ). 


respect to r change with time without changing their form 


they remain similar to themselves. This type of motion 


in which the distribution of the flow variables remain 


similar to themselves, (i.e. similarity in the motion itself) 
with time and vary only as a result of changes in scale is 


The function p ( r, t ) can be written in the form p(r, t) 



» ^ 
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- tt (t) P (r/R) , where TT (t) is the scale of the pressure 
and R(t) is the length scale both depend on time in some 
manner, and the dimensionless ratio = P (r/R) is a 

function of new dimensionless coordinate *1 = r/R. Multi- 
plying the variables P and by the scale functions Tr (t) 
and R(t), we can obtain from the function PC'] ), indepen- 
dent of time the true pressure distribution. The other 
flow variables, density and velocity, are expressed simi- 
larly £15] ^ 

It is a natural question that what requirements must 
be satisfied by the conditions of a problem in order that 
the motion be self- similar. 

Dimensional analysis is used to answer this question. 
Since the dimensions of pressure and density contain the 
unit of mass, at least one of the parameters in the problem 
must also contain a unit of mass. In many cases this is 
the constant initial density of the gas -f^, which has the 
dimension M L ^ • Let the parameter containing the unit 
of mass is a. It can be assumed that its dimensions are 
aj M L T . The dimensions of the functions p, f and u 
are, [p] = M L ^ , [-fj = m L"^, and [u] = LT“^ , 

we can without any loss of generality, represent them in 
the form suggested by Sedov [ # 



k+1 . s+2 
r t 


^k+3 . s 

r t 


Where P, G and V are 


For self similar motion it is possible to reduce a 
system of partial differential equations to a system of ordi- 
nary differential equations for new reduced functions of the 
similarity variable , R - R(t), The boundary and 

initial conditions of the problem are made dimensionless and 
in term transformed into conditions on the new unknown func- 
tion of '1 • This simplifies the problem greatly from the 

mathematical standpoint and in a number of cases makes it 
possible to find exact analytic solutions. 

The problem of a strong explosion represents a typical 
example of a self— similar motion. This problem was formulated 
and solved by Sedov and succeeded in finding an exact 

analytic solution, to the equations of self-similar motion. 

The same problem was also considered by Stanyukovich in his 
dissertation (jlt) and by Taylor , both of whom foirmulated 

hhe equations for the problem and obtained numerical and not 
analytic solutions. 

The parameters in the problem of a strong explosion are 
the initial density of the gas fo'^ and the energy of 

explosion £ I ^ The energy E is always equal to the 

total energy of the moving gas, and as a result an energy 
integral appears in the problem* These two parameters can not 
be combined to yield scaler with dimensions of either length 
or time. Hence the motion will be self-similar, that is, will 
be a function of a particular combination of the coordinate ' 
r (distance from the center of the explosion) and the time t. 
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The initial pressure and speed of sound p and c in the problem 
of a strong explosion are assumed to be equal to zero, and hence 
these quantities are not parameters of the problem. So the 
quantity r/t can not serve as the similarity variable. In this 
case the only dimensional combination which contains only length 
and time is the ratio of E to fo , with the dimension 



[E/*J 


_ ^ ^2 ^ Hence the dimensionless quantity 

/ fo 


/ To i 

1 E1tV 


can serve as the similarity variable. The distribution 
of pressure, density, and gas velocity can be expressed as 
functions of. one dimensionless variable ^ » 

s 

8. CONCEPT OF SELF-GRAVITATION 

A fluid can be referred as self-gravitating, when the 
fluid mass included is large and isolated and the gravitational 
attraction of other parts of the fluid provides the volume force 
on any particular fluid element; for example, a gaseous star. 

In the case of spherical symmetry the effect of all masses 
on a particle (point mass) at a distance r from the centre of 
symmetiry is equal to the force of attraction by the point mass 
placed at the centre of symmetry and having a mass m (r, t). 

The behaviour of gravitationally interacting gaseous 
masses forming a star is given by the appropriate equations of 

motion([l2], [lej) . 

The continuity equation for radial gas motion with sphe-’ 
rlcal syimnetry is. 

It 4- u 11 +- f { -fx 
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I i 

For the invlscid gas, momentxam equation taking gravi- 
tational forces into account (according to Newton* s law of 
gravitation gravity g at r is ) can be written as, I 

i 

alt ^ 3^ -f TS ^ ^ 

'■ ' ? 

Where G is the gravitational constant and m (r, t) is I 

J 

the mass of gas with in the sphere of radius r, j 

• I 

To determine m we use the equation, | 

which expresses the fact that the mass of a spherical shell 
of radius r and thickness is 3A. 

and the energy equation for adiabatic gas motion. 

These are the four equations with four unknowns p, 

-f , v and m, 

A very large no. of papers have been published by now 
in which analogous self-similar solution were abtained to 
explain the adiabatic unsteady flow in self-gravitating gas 
and analysed for systems of partial differential equations 
encountered in various problems of Astrophysics such as 
internal motion in stars, motion of nebulae etc. 
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CHAPTER-II 

A SELF-SIMILAR FLOW OF SELF-GRAVITATING GAS BEHIND 
A SPHERICAL SHOCK WAVE IN MAGNSTOGASDYNAMICS . 

1. INTRODUCTION! 



The explanation and analysis of the internal motion of 
stars, the theory of explosion of novae and supernovae stars, 
the motion of gaseous masses at high relative velocities with 
shock waves in a gravitational field are the problems which 
are of the great interest now-a-days. Carrus[l9] and Sedov[l^ 
obtained the numerical solution for the self-similar adiabatic 
flow in a self-gravitating gas. Starting with these results, 
Ryazanov |[2o] obtained a particular analytic solution. But this 
solution does not describe the flow behaviour in general. 

Singh j^ 2 l] bas discussed the self-similar adiabatic flow of 
self-gravitating gas in ordinary gas dynamics and has obtained 
numerical solutions. 

In this paper, the propagation of a spherical shock wave, 
in self-gravitating gas, in a non-uniform atmosphere taking 
account of magnetic field effects on the physical parameters 
is considered. It is assumed that magnetic field effected 
only a portion of sphere enclosing the origin i.e. the point 
of explosion. Similarity principle has been used to reduce 
the equations governing the flow (inviscid) to ordinary diff- ■ 
erential equations and it is assximed that the density varies 
as the inverse power of distance from the point of explosion. 



The total energy of the flow increases with time. 
equations of the prdrt.-rm 


The equations governing the adiabatic flow and field 


are. 
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( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


The flow variables just ahead of the shock, denoted by 
suffix zero, i.e. , -f^ , he, , are. 


\x^zz a 


-K 

> -f. = ,■ Vn, = 


( 2 .^) 


Where R is the shock radius. 
From (2.5) 

■1^ = 

at r = R ( at shock ) 

- . ^ - '^ -2- T=. 4 -TT f e, <<, = M TT 


as -fo 




so, ahead the shock. 
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(after integration) 
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Where 


) i- ^ 


and A, c, "<v j (S 


are constants 


The Rankine— Hngoniot boundary conditions of the shock 
surface, following Whitham [22] 
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Where, 








( 2 . 12 ) 


V * dR/dt is the shock velocity and is given by the 
ratic equation. 
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M ch number M and Alfven* s Mach number are given by. 


>' f. 

H =: — ^ 

iK 




K 




( 2 . 14 ) 


The permeability constant of the medium {>>- ) is taken 
to be unity, 

3 > SIMILARITY TRANSFORMATIQWS 


For finding out the solutions/ some similarity assump- 
tions are made by writing the unknown in the following form. 
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UhJ , -f = 
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( 2 . 15 ) 


Where 


H = A 


b 

i = A % ■ 


Ix+aJ/i. 

a ^ 


nC'ij 


( 2 . 16 ) 


and K, A , a and b are constants and are to be determined 
from the conditions of the problem. 

The total energy E (which increases with time) inside 
the shock wave of radius R is given by 


£ = H tt 
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i - P- 1 - = a*!- M , 

^ k J ^ ( 2 . 17 ) 






Where B and q are constants and r^ is the coordinate of 
inner expanding surface. 

In terms of variable vj , the total energy E can be 
expressed '-as. 


Sis 







^- 2 . 2 . A- 2 - ^ i 

it iig^+- A. % Tpp ■’■A t M ~ 


XKt-'SL j|.x v 'i. tjL 

- 6 i/i vjil j A. 


K+H A-l / 5^ . JL.4.M 

4 % ±0 y^, 


iKt-s iA r ^ 

- GiA. ;t Wii. k 4/1 l 


Since# 


also# 


-i- -^- di 
a. ^ » ' 

j<±C -(K-i-rj'o/(x. 


ait+iT - k K-f- si *>/ 0 ^ 

A = >1 t 


Substituting these values we get# 


_ MTT 
~ . a- 


K+S" A-2- Xf^ . -t-H ) - 

A- t / 0 it •t Y-t y 

vl 




where 'l'*' and vj being the values of »{ at the expanding 
surface and shock front# respectively. 


- G^( 4 A *1 t 


Hence by (2. 17) # 


. _ . ' 



We have 


constant at the shock surface 


So velocity of shock is# 


(ahead the shock) law is 
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From equation (2,19)# we assume# without any loss of generality 


Using (A) 


4o SOLUTIONS OF THE PROBLEM 


The equation (2,16) is used to reduce the equations 
(2.1) to (2. 5) into the following forms# 





and/ 


Since 


as, _il. 
at 


,2. . J<.+-2 A~l 

^ ■— A % U(iij 


= t'^ u(^tj -t- 




K+ 2 . A-l 


-^'(lj Ul'j 


I 

«-lA u(yji(y 


using these values equation (2.1) reduces to. 


^1. a Wi j u (^ 


OL'j 


OL'^ 



\ W>^a'ClJ+ ♦ a^.nt'ljUhj t 







(^K+ijA t X rcy ^ 


IC+i A-Z 


where 


^ =1 (K+^j A ^ + o-i It'^'^Vcij 

H - M(^j 

= i^) '^Cij +- M a'=Z= t'Z*-' 

Just ahead of the shock we have, 

2 -r[}/^^ 6 i a-a.K ^ - 2 S 

I t- ^ = “< 


We assume that motion of the shock wave satisfy the 


power law 



■fo = AA 


A = 


and magnetic field distribution law is 

Kf, = C ^ ^ , 




2 ^ 2 . 2 A 

nf-(t, ^ <h ^ 


(t-#;R 


, , , , , . ; %’ t '; •; ^ '. / 


V 












» A 

r4Ni L 4- JbL- 

a 


^ -X. r%* f 


a.TT ft- 


substituting the value of a,b,K, and X 




Jii + 

■■■JC . 


(?- tjC'^-'j too ^ _i_ r -L- - 

vtTft ~ o.s'K®- ^C"<-') 


- -1 4-J 


Where L- [ 1 ^'^ 


J__ - 


Equation (2,3) is. 


Since, 







,*Va 


oiij -» “'i i' ^’(ij 


using these relations equation (2,3) can be written as 
Z?"^' NJ^J -h *=1 ^ ra’tij +■ 

JSu 4- 1 ^ JE, f-i J — 2 

-»- NtiJ ^Cij + u’Cij ^>(.>1j •< 

tl JL - z . , „ A 

+■ ^ it’- ° 

or : 

-4- + A*|u*(iJ Ufij 0 

substituting the value of a, b, A and k we get, 

-M ^.' 1 Kj' + (l- ^ MU - riw'u + aMU -S-lu'-u 




)l 

ri -ij- 


-j- au - t- (b-ij 


-U =. ti 


Equation (2.4) is. 







substituting the value of IX -’T equation ( 2 . 5 ), 



' Ki* 2 l 1 !✓ 4 . !> i 4 

(l<+a) /t loCij -t cx..j A ^ ^ (ij 


or 

(K+3j w(ij +■ 0 ,*] w'dj = Mir Atij 



or 

tK4-3j + ^ M TT 



put 

f< - ^ K and • ^ ^ ^0 g^-t^ 




S'l -S" " (S-tJtM-n -ij- =0 


(2.26) 

The 

shock conditions (2.7) to (2.11) are transformed. 



^1 — "^X 
f u(,J = -1^, 

(by 2.7) 



when v| - , R = r 


■ , '! ’ '' 


X ^ “X t 


* „ ‘ i ’'J 


un.j = - K - 


" ■ ' ' 

or 





u(i.j= 


( 2 . 2 tj';' 

By 

( 2 . 8 ) , 


aiii^iaii 

iKiiilMl 


-f, = 






K; 
















substituting 




1^ = 


By shock condition (2.11) 


A,UC^J = 


(2-^; 




A _ 

t2-<j 




mCIo) 

A 


In the form of non-dimensional equation (2.22) to (2 
to be transformed. 


By (2.22) 


-o 


or 





and prime denotes the differentiation with respect to ^ 


The jump conditions at the shock are 


Equations of motion are. 
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5. NUMERICAL RESULTS FOR THE SIMILARITY SOLUTIONS OF THE PROBLEM 


have been used in the following non-dimensional forms. 


K+-t 









(A 


Numerical integration was performed on a computer 
using the Runge-Kutta prograiwise for the three eases *- 
V » 2 and <>( = 2.5 also for the simplicity it was assuiw 
: ‘I - ‘1„ :=. I • The other constants 'are 





pressure and density, temperature 
shown In the flaures No. n ^ _ r 


6* CONCLUSION 


In the problem of propagation of a spherical magnetogas 
dynamic shock wave, the velocity and mass are maximum at the 
shock surface and decrease on moving in opposite direction to 
the point of explosion from the shock surface whereas the 
density, pressure and magnetic field are minimum at the shock 
surface and increase on moving in opposite direction to the 
point of explosion from the shock surface. 
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CHAPTER^III 

A SELF- SIMILAR FLOW BEHIND AN EXPONENTIAL SHOCK 
WITH RADIATIVE HEAT FLUX 


i. INTRODUCTION 


The problem of the self-similar propagations of a 
shock wave in a exponential atmosphere taking into account 
radiant heat exchange has been discussed by many authors 
viz. Raizer |~23 J • Zel' dovich 24] « Hayes [ 25 ] etc. It 
is the exponential character of the medium (density have 
he exponential law behaviour) which ir.p .rts to the motion 
its special features. ** 

Gusev [25 ] Ranga Rao and Ramana[27] have studied 

the problem of unsteady self-similar motion of a gas dis- . ; ■ 
placed by a piston according to an exponential law. Verma 

and Singh ^ 28] Singh and Shrivastava ^29] have consi- 

■ 

dered the problems of spherical shock waves in an exponen- ,, 

■ ■ ■ ■ ■■ ' ■ ■ ■ ■■■ ■■■■ 

tially increasing medium under the lew uniform pressure 

including the effects of magnetic field and thermal radia- 

.. .. ^ 

tion# respectively. ; . 

In this paper, a self-similar flow of a perfect invit-',;,';' 
cid gas behind a spherical shock propagating into an uniform 


Im 




exponential medium (which is assumed to be at rest) is 
considered# The surface of contect discontinuity movesS’:' 
time according to an exponential law. the similarity solu- 
tions have been obtained taking radiative flux into 






deratione. Radiation prasanre#' radiant energy#, sol^r grav 


■ 










2« EQUATIONS OF MOTION AND BOUNDARY CONDITIONS 


For inviscid perfect gas, the differential equations 
governing the one-dimensional unsteady flow with radiation 
flux, behind a spherical shock are. 


where being the gas constant; and F denotes the heat rxwc 
Assuming local themtiodynaniic equilibrium, and taking 
anrt* «? diffusion anoroximation, which is. 


where is the Stefan-Boltaman constant, 


path of radiation# 

Following Wang ^3oJ 
erature and density, is 


which is a function 




moves with time accor- 


^ (^>o) C3.7) 

and since the flow is self-similar, the shock will also move 
with time according to an exponential law 

a - § ( 3 ^ 3 ) 

where r is the radius of inner expanding surface, R is the 
shock radius. A, b, and m are dimensional constants. 

The shock boundary conditions are. 


miiitsie aujjscxipcs u ana i aenote the region just ahead , 
and behind the shock, respectively, and v is the fluid vel<i< 


3. SIMILARITY TRANSFORMATIONS 


The similarity transformations for the problem under 
consideration are, i ■' « 






U^ = U-UllJ > 
■f = J 


80 

( 3 . 14 ) 

( 3 . 15 ) 

( 3 . 16 ) 

( 3 . 17 ) 
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using (3,15), (3.16) and (3,17) wa get 




c/to .K a^-t-S <^-(1 0 j f p'^ti 

3PW+P h j T'ty 3 ^ H 




- ^ g-VP'jHO .>V-MyVa' |jj 




using (3,18) and (3,19) 

3 -p M+ (i 




- M 


f^-^l 


{3-. 20) 


with P~- 2 .. f oi remaining arbitrary, {o^<=< ■5-a- ) and 


nJ = 


4o->M c Mo 'f<' 

3 -p H t- p. 


— a dimensionless radiation parameter. 


4o SOLUTIONS OF EQUATIONS OF MOTION 


The equations (3,13). to (3.19) are used to reduce ill 
the equations (3,1) to (3,3) and (3,20) in a new form- 


f-= 

Jl- = -fo Tt 




■fo ^ Clj ' ■* 


JU 





also 


so. 


and. 


il 


li 


Cij 




u. 


il 




u. Jl 

^4. 


IC 

a A 

h ati. 

^ aA 
2.i^ 


VA1 Gi^CnJ v/tij 




a 


il 

iA 




= ill- 6*tij '/tij 

Substituting these values in (3.1) 


ii_ ^ ^il .+. 


!>-k 


aA 


+ f 






or. 


or 


-t Gi’v) -t ■^- 

*tv‘j 

Gn =. - 

lit / .M kii 


o 


( 3 , 
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Equation (3.2) is, 

UL - 


"at 


since/ 


and 


and. 


au 

it 


a*. f 3 ^ 

a =: d v/tl j , 




^ (it 




a^^ 

at 


R 


^*1 _ _ 'mA 

d 


30 . 

at 


a/i- 

3 U. 


~ _ VA^ «! 


u. 


a/i 


O'*" ,V 

X '/h; 


a — vsitJ- v'(,^j U(tjj 




Jt- 

»A 




J- at» 
t 3 A. 


-t- llL 
■f aA- 










ycij 

Oltij 


(G^ t-Kj 



o 


Substituting these values in (3,2) 


-+• iJ- Vt*! j +■ v^i> y(.jj ■+ I,v\v> 


or 


or 


- vi v' +- V -H \/*V/ +- JE_ _ ^ 

G| ' 


J&L 

Cr\ 


•iw - M - V V 
&i [ v’(i-vj - vj 


( 3 . 22 ) 


Equation (3,3) is. 


it 


u 


j/V 


-it - -XiiiL ±L -H y-' 

^ at f T\ 




Hij - 

— ■(. [ 'P(i( - 

[’-'PhJ - .)T'(„] 

slibstituting the value of P* from (3.22), 


> 1 ° 

d-t 


- 'l&i iG,vJ 


1 


j±_ =- 

aA. e 






'f ts 


u [^Gj (i-vjv' ~ CHVj ', , 


’■i' ! 

*’ ’ ' ' f' I 
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and. 


J±- Ji_ 
f ai- 


- Hfij 


substituting the value of from (3.21) 


so. 




a -f 

i ifx. ' 


i^vdl 




5 ^ T^v/ (iv'-f 


“li l-'^j 


aA 


UVj 


aA 




fo l> ^ ^ ^ (ij -h Q*(lj 


= (aa -t 


illLL JL. 
A*- 


(a"f; 


. -I. ' f' 
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61^ - f - V - ^^^'•^11. ] 

L 'ili-''; J 


or 


AO.""*' ^ _ «! - _ ^''■/o.JX 


»J 


lU-y) iti-'^j 


or 


_ - ('?/6fji/ - a.(p/cju 


or 


V' [ ~ v^ Ci-vj - {‘i-'/j 


{i-sj 




ft Gj 


0‘</ 


►/'P 


''] 


or 


(i-vj 


a. V 


. . r a»Gi ' ..1 

('W(h) - h ( i -V [-TTf” " y 


1 [ ti-'^J’-- P/&»j 


(3.24) 


• • j 
' ') 


The shock conditions (3,9) to (3,12) are transfonmed into the 
following forms, ( at = 'J, ) 

Boundary condition (3.9) is. 


flpP 

; ! ,i 


u, 




: ; .‘it/' 


v> VI tl) = [ » - fH®-] 


' '« - •' 






iS* *j'’ 

Condition (3,11) is, 


Thus the transformed equations of motion are, 


^ f ( I “ -^)J - ^i^ii*-y\^J - = lY‘tj 

UV'fJ 

and 

^1 _ ay-p/Gt - iCi-'/j [^*’' *^^/m7» - ''J 
M [I'l-'^J^ - P/6,J 

And appropriate transfonned shock conditions are^ 

V(tj 

'Pf.j 

and 

61t>| 

where prime denotes differentiation with respect to 

('•i 

5. NUMERICAL RESULTS FOR THE SOLUTIONS ■ - 

^ — — , 

' ,4 

For numerical calculation, it is convenient to writ.#)^, 
the flow variables in the following non-dimensional forms,'; 
(using (3.14) to (3,17) and (3.25) to (3.28) ) 


t ' . 


YM 

I 




. j 

i'i 


-i- 


If 


G| 

^OJ 




(3.30) 


JE- 




T>(i; 


(3.31) 


and 




4-r-4 




■]t. 




f^. 


au; 


(3.32) 


The numerical integration was carried out until the 
kinematic condition is satisfied and was performed 

on a computer# using the Runge-Kutta programme for N = .10 
and N = 100. The other constants are .M » 20, y=-»*4 and 
I . The nature of the flow variables is given in 
tables I and II for N “ 10 and 100, respectively. 


f 

. 'tM 

f/Vi''' 


6. CONCLUSION 


I.n 'the problem of the propagatioii of an exponeipitiaX,!J' 






.shock wave with radiative heat flux# there is very minute 


mm 


variation In velocity throughout the meaium, while the. . 


variations in pressure, density and heat flux are ireiry 






minute in the beginning (near the shock surface) and then 
rapid in the last (apart from the shock surf ace) # For bvth 
the cases N * 10 end 100 the velocity, density and pressure ' 



air© rEiinirnvitn at the shoc)c stitface ( i»e« at '1 ~ i ) aRd increase 
on moving toward the point of explosion from the shock surface 
whereas the value of heat flux is maxirnian at shock surface 
and decreases on moving toward the point of explosion from 
the shock surface. . 
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J. H. iilNtiM AN« H. f, VtttlWAKAKMA 


TAULEJ 
A « 10 


f.UUUU 
1,02160 
1.04273 
1.06344 
I.0H377 
1.1037K 
1.I23S2 
t I4.i0.l 
i.l623S 
t.U15l 
IJJ0J2 
1.21042 
1.2.1KIU 
1.236117 
1.27342 
1.20200 
1.31227 
1.33034 
I.3467U 
1.36673 
1.3K46lt 
1.40240 
1.42017 
1.43771 
1.43311 
1.47234 
1.44040 
1.306244 
1.52206 
1.53042 
1.SS56S 
1.57161 
1.54720 
1.60263 
1.61766 
1.63227 
1.64643 
1.66013 
1 67.1241 
I 64360 
1.607340 
1.70416 
1.71746 
1.72623 
1,73204 
1.737444 
I.730U01 
1,7.3664 
1 72644 
1,70351 
1.671 16 


l.OOUUU 
1.01021 
1.04060 
1.06421 
1.00014 
1.11444 
1.14034 
‘ 1.14244 

H. l »|,j3027 

1.23472 
f;*' I ■201 30 
1J4740 
^ I..10K24 

14.320.1 
1.31242 
1.57725 
1.64700 
2.72537 
1.41010 
1.00220 
2.00607 
2.11035 
• 2.24324 

2.24324 
134133 
171676 
1014103 
2.13746 
3.30161 
3.62451 
4,01763 
4.40646 
4.461 14 
5.30767 
6.036042 
6,40710 
7.747104 
4.21040 
10.174V0 
12.23600 
14.72440 
14.06007 
22.73064 
20.556 45- 
40.07473 
57.33726 
40.43122 
150,0425 
361.0603 

I, 467.0463 
0727.42U4 


13X1000 
0.07714 
0.06630 
0.06720 
0.07060 
1.0O2O4 
13)3334 
1O7700 
1.12014 
. 1.14004 

1.23752 
1.334 71 
, 14200} 

‘ ' 1.31676 
1.62241 
, 1.74023 

* U700$ 

1 2.01267 

117240 
134046 
154601 
176524 
2311063 
204612 
’ 13M33 

3.04771 
4,34674 
4210224 
5.32302 
5.02425 
, 6A2467 

744743 
441170 
0.53677 
1002013 
1130102 
14.627^- 
17.13332 
20 ,14331 
2442312 
20.40007 
37.01473 
473)1767 
6143244 
41,23440 
114.3664 
141.2306 
m3204 
64424)04 
216*44617 
10210.132 


1310000 
0.374674 
010146 
-016242 
-044417 
-07wo7'il 
-13)7671 
- 1.33277 
- 1.62344 
-1.00247 
-110017 
2 40403 
•242123 
-1.14113 
-2.570023 
-4.021076 
-4 320 14 
-3,II4«03 
-172301 
-&46024 
-7.30*30 
-aiTOSO 
-041077 
, ,*c»»030 

•’'■"S'- 1127438 
-143IW03 
-t6.1«tlO 
-UJ(1I74 
&213I4443 
>-23.30000 
-30.11200 
-33JM343 
-4144363 
-314*242^ 
•*6241*34 
■;^7*i3H)2 
4#‘«03i7737 
T. ' -H7.2I34I 

■ '■ * 144.06023 
#■- -I013I804 

-341473SS 

^ 6464143 . . 

■ *tU42442| 

-36147271 - 

■■'■ -2«»7.7*0| 
-6t26J032 
■'>|6mi04 ■;' 

-638 ia »7 


A !!itl.»-S)MU.AK l-tow llUilNU AN tmiUMJIUjU SHIA'K. i 

tabuu 

Nsm jr, 


131 


I uu 

OWi 

U‘A< 

uvj 

{I'M 

II Hh 
. II H(i 

U.M 

ua3 

IIIUI 

ii7lt 

U.7o 

g.74 

«.72 

UTU 

U.W 

UMi 

U.b4 

U.(i2 

UbU 

gin 

0 ill 

UM 

US’ 

UXI 

U4« 

U.46 

(M4 

(MS 

IMil 

U.Ji( 

{i.M> 

OM 

UJ2 

U Ml 

U4M 

o,2t. 

(I .’4 
U 2 .' 
lijO 
uU 

U id 
U U 
Ul.' 
Will 
ilWh 
u IA> 
(Ilk 
uu.' 
U.UI7 



I.UUUUO 
I.U2UH4 
1.04161 
1.062.12 
1.0K2y7 
1.10.i5.i 
1.12407 
I 144.^4 
I I64S^ 

I l»S29 
1.20iJV 
. 1.225113 
1.24601 
1.266 IS 
1.211623 
1.30625 
1.32623 
1.34616 
1.36603 
1.36566 
1.40563 
1.42535 
1.44502 
1.46464 
1.46420 
1.50371 
1.52317 
1.54256 
1.56190 
1.56! 17 
1-60036 
1.6IVSI 
1.63657 
1.65754 
1.6743 
I6V522 
1 7I.IW 
1 7.1246 
I 7M)I(',I 
1 76'! 17 
>1.711726 
I 60516 
1.62265 
1 64023 
I 65726 
I 67 163 
I 66V6I 
I SKWy3 
I y|670 
I y2y6o 
I yiii2 


1.00000 
1.01097 
l.tM171 
1.06533 
1.09095 
1.11671 
1.14676 
I 16133 
1.21656 
1 24571 
l.2V<iU4 
1.34065 
U6946 
1-44226 
1.49974 
1.56233 
1.63065 
I.7053S 
1.76722 
1.67714 
1.97616 
2.(61555 
2.20673 
2.34143 
2.491 72 
2.66006 
2.64961 
3316367 
130706 
3J6529 
3.90535 
4.27613 
4.70905 
5.21692 
5.62530 
65,'M45 
7.44224 . 

6 5.1677 
9.91 >73 

11.77116 
13.9660 
17.0799 
21.4299 
27.7259 
37. .192 3 
53.3690 
62.6622 
145.632 •, 
325.764 
1295.614 
10565.264 


13)0000. 
0.99726 
0.9963$ 
0.99726 
100017 
1.00514 
1012)1 
102164 
103391 
10467m 
1.06643 
1.06735 
1.11173 
1.13969 
1.17217 
1.2U699 
1.25079 
1.29606 
1.35147 
1.41164 
1.47936 
1.55562 
1.64144 
1.71613 
1314719 
1.97045 
2.11007 
2.26670 
X44956 
2.65660 
2319475 
3.17019 
349075 
3.66650 
4 310SS 
414021' 
5.47679 
6 25623 
.7 22333 
6.43656 
9.99963 
12.0536 
14313662 
IK7498 
24.5193 
31.5657 
49.175 7 
79.9317 ■ 

156.921:;'. ■ 

465.IS1 
2577 71H 


1.00000 
067336 
0.33063 
-0.02956 
"0.40831 
"0.60656 
1.12321 
-1.66149 
-2.16413 
- 2.669 16 
-3.31444 
-3.79946 
-4.42914 
-5.109U 
-5314617 
-(iA47W 
-7.52315 
-6.41276 
-9J3913 
-107070 
-12.0040 
-13.4511 
-15.0735 
-16.9017 
-16.9727 
-203)6 
-i4m4 
-27.1J0I 
-3ft7W 
-3451907 
-39.9654 
-45.6706 
- 52.9440 
-61.4993 




' . :.r 


~ 101,1491 • 

121645 
")466«9 
-I64.2J7 
•212,471" . 

•>«7439.,v 

'-$45,143 

- 761,996 ' 

• IW.IlS 
-3769,76* 

. “41556.363 " ' 

37945104 
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A SELF-SIMILAR FLOW OF SELF-GRAVITATING GAS BEHIND 
A SPHERICAL SHOCK WAVE IN MAGNETOGASDYNAMICS 

O. VERMA* andV. K. 

DejfortmerU of Mathematics, University of Gorakhpur, India 


(Received 7 Octoter, 1985) 

Abstract This paper considers a spherical shock, in a conducting gas, of self-gravitating gas propagating 
in a non-uniform atmosphere at rest. Similarity principle has been used to reduce the equations ^venting 
the flow to ordinary differential equations under the assumption that the density varies as an inverse-power 
of distance from the point or explosion. The total energy df the wave is variable. 






IiitrodiH:ti0o 


The unsteady motion of a large mass of gas followed by sudcteaa release of oiergy results 
in novae and supemovae. The explanation and analysis of the internal motion m stars 
is one of the basic problems in astrophysics. Nummcal solutions for self-similar 
adiabatic flow in self-gravitating gas were obtained by Carrus et al (1951) and Sedov 
(1959). Starting with these results, Ryazanov (1959) obtained a particular analytic 
solution. But this solution does not describe the flow behaviour in general. Singh (1982) 
has discussed the self-similar adiabatic flow of self-gravitating gas in ordinary gas 
dynamics and has obtained numerical solutions. 

In this paper we show that the magnetic field has a significant eflect on the physical 
parameters when the gas is self-gravitating. The motion of the^ shock wave is assumed • 
to satisfy the power law 

^ ai) 

where A and a are constants and R is the shock radius. The total energy of the flow 
increases with time because of the pressure exerted on the gas by an expanding surface* 
Therefore,’ ■ 

£=s/v( 9 >o), ..viv'ai) 


^ ( 1 . 2 ) 

where £ and 9 are constants and E is the total The maga^ Md dishii»Akm ’ 

law is 

ho^CR-^ (P^O), , (1.3) 

4 „ i , 

where C and p are constants, and the values 0^9 and j? ate to be d«*«i*w«l laKr, Tike 
magnetic fieid is directed tangraitial to the advancing shook ISrcwat. ‘ ‘ r 

f i; ’ ^ ■' : 

Professor Verma passed away on 4 April, 19S5. ^ v?-'. 1 1 r| 

** Supported Hy CSIR, New Delhi under grant No, 7/57/M7i^t./EWR4. _ ^ ^ 
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The flow variable just ahead of the shock are 


Mo = 0, mo = 


— 2j^ y 

" («-l)(3-(i) 2^' ■ 


where 


1 4* ^ = a 


and viscosity is n^lected. ’ '* ' ■ 

We investigate three types of the first having total energy of the explosion to 

be constant (a = 23), sk^ohd having constant velocity of propagation of shock 
waves (d « 2) and the third having neither constant total energy nor constant velocity 
propagation of shock waves (ot » 1.5), An idealized magnetic field is considered for only 
a portion of sphere enclosing the origin - i.e., the point of explosion. 


2. EquatioiB of the Frobleiii 

The basic differential equations governing the adiabatic flow in a seif-gravitating gas are 
given by 

dp 1 d \ _ 

- P‘) 

du du I dp H dH Gm „ , 

— + M — + - — + - x — + — + — = 0 , (2.2) 

dt dr p dr p dr pr r^ 


dt\p 


dH dH „du HU „ 

— + u — + H ~ + = 0, 

dt dr dr r 




■ S::dr 


4npr^ , 





where r, 4 m, «, p, p, and y are radial distance from the centre, time, mass v 

in a sphere of. radius y, velocity, density, pressure, an^ ra«o of two 11^1^1,1 ‘,!,i 

respecA^^tod ^ G represents the gravitational constant. . . ^ '• ' ' i” 
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( 2 . 15 ) 



where ij = r“f*; while K, A, a, and b are constantse and are to be drtermined from the 
conditions of the problem. The total energy £ inside die shock wave of radius R is given 


In terms of the variable r\, we can express the total cnwgy as 


- GiVil ■ r " 'J d(J »• af ; 

»}(, and ij'^ being the values of r\ at the shock front and cxpandiog swface. 
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Equation (2.17) yields 


This choice fixes the velocity 




•t*'2.00 


Fig. L VeWly dfetiiltifoa 
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where ii is assumed that initially /? * 0. From BqpiatiOB ( 121 ) it can be seen that the 
value q = 1 corresponds to uniform expansion the sjrfiwe. , ? 

•i|‘' >■ "'I V' ^ ^ 

3. Sirfutioiis ■ 

The condition inside the wave will be obtained frenn Equations (2.1)-(2.5). Using 
Equations (2.14), and (2.20), the equations of motiem are transformed to the 

following non-dimensional forms 






, "V'- 'V ' 
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where 


at the shock 


0.90 




Fig, 5. Urns 




are given by 




0 'i ^ , 

Results and Discussrai > 


Similarity solutions of the p-olrlem of 


spherical mt^etogasdynamic 


p _ P(ri) 

Pi nJ Pi>h) ’ 

^ V%/ Nino) 


(4.3) 

(4.4) 



Numerical integration was paformed on a DEC system 1090 computm^ testalied at 
I.I.T., Kanpur using the well-known RKGS programme fOf li» a - 1.5, 

a = 2, and (X = 2.5. The other constants are f «• «Bd a • ~ 5. 

The variations of vdocity, density, pressure mapdic and mass have been 

illustrated graphically, sec Figures 1-5. 

At the shock surface, the velocity and the mass are mwtiianm, 
move away from the shock surface, see F^urcs 1 ad4 S, !%*» 
that the density, pressure, and mt^etic fidd are 

, . 

increase as we ni.ove away from the shock smlliee^ ' 1 


’ i| '1^ ' 

^ Pf?^ WmfP '^rw^pppi' 
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SELF-SIMILAR FLOW BEHIND AN EXPONENTIAL SHOCK 

^ W»TH radiative MEXTFmXi U . 4 


' / “ (.(SrtSfw I- ^ ^ 

^ ^ ^ J. A SI N^H Kiij^ P. R. VlS,»iyAKAjRMA 

Uoi'AikCMk < Wfcw, iM-ki. 


tki, a.U.Rt f’uili/fiutumi/ Cutkyi^, Uvhhi, Juniipuf, India 

h-;-- 

' -'f;' ■ 

(kcLcivcU 2 FJrtiniary. W3) 


A^viLFatl. Si‘f( SUJUj4/ Tof UIIif’‘iillUCUSii>HUl UU»t;‘a(Jiy u| a JH^tllvCt g4;v fpclUHti I4M CAliUllCiJli^J 

khiH,K Hito 4ID itm uik| i|*i4 urn vuDAnitfii- 

liwii. 'I hc u;i4l energy of Uic Wilve k vg|iul4‘, * 

I*' ‘ N"4,' i>» / 1 / ' 1 '*'4'fe - ‘ * , ' , . 'M‘ 


' I ,1 ' * - 


/ ", ;' ■;, <' ' :' y, ■;'; ■ ■ ’ • I. lBlfa4ucti«« ' 


* ' 4 ,, ‘ , 


Steady wlf-simiJar motion of u pi» Uivpiaccd by u piston uccorUitig lo tin exponential 
law. Vcrnia and Sinijli f 1979/ «nd Singh and Srivasiava (l9S2;j have cutWdered the 
[)rul)ldins of spherical sliock waves in an exiK>ivc»iially increasing tnediuni inalcr tlw 
low uinliHiii tticssiirn iTiduding tlw effects of iiwvneiic field and ‘thciniaf widittiom 
Test>ecli'vely. 

I he prcseiK paper deuJa with u setf-siniilar lluw of u tseriix't gas behind an expo- 
nential shock driven out by an expanding surface, i.e., the surfaw of oomaot dt»- 
(.oiittiHiity which moves with time uevording to an extHmeitliul law, Tlie vhotik waves 
piopagale in an uniforin aluiuspheic which rs astiumed to be at rest. TTjb tioutafity 
solutions have been developed when the radiation heat tlux is more intixirtant tfias 
the radiation pressure and radiation energy. Ihc oiliet assumptions include Uw |as 
to be grey and opac|ue and the shock to Ik- trauspariitit-aod isitdiefoial The total 
energy of the wave vanes as the cube of the shock radius. Viscosity and tlw solar 
gravity have been neglected. 


A-i% 


2. Equations of Motion and Ikauntarv Couditions . , 

■ • 

How behind a spherical shock a»x »' 


The equations of How behind a spherical shock a>x 


dt •"P'5; 


(ur) s» 0, 


I ’ ' * Vi'' 

f, , . (H 


rfu I ap ^ 

it _ If it 'i fcil i. 

dt /) fit ~ f** df 

f-w , , 






.,7" 

'■'Vr . vf 






A xtn^pkyMa mti Siimv c» iai4i4, m4^ : l. 

i hy fii*M MIMm ^V.s IJwdmk md ^ % 







uppruxiiiiutioit, wc iiuvo 



wliuic fcj^ IK llic Stafun llull/mniMi aiiikluiit^ e. lliu vuliK'ily u( liuhl. iv liiu iiitfuii fiwt 
pulli lit luUiuliiut, Ik a luaulluit ul dvpkUy uuit ii;u»|wiutuie. t'ultuwiog Wu(i)|(i06t!^ 
wc.tukc 


Pu . «, unit /f UiiUK aiiuiluulii. 

1 tic innpr expanding iurfucc muva )yilh time acvufding to an exponential taw. 

f»/lexp(m<), («i>0) ^ t?) 

und kince wc have uiikumcd kclf-kimiluriiy, the sslUKk witi abo ntuve with time ocoird* 
mg to an exponential law. 

KmBexptwt^, Cl) 

whereT in the rudiuN of inner expanding kurruev, R in tiw nhuck radtut, A, and m 
ore dimenxtonal constants. ■ 

' The disturbance is headed by an isuthcftnul shuck; ihererure,tlw boundary con* 
ditions are . 


where subsaipts 0 and i denote the t^gNa and behimt the 

shock front, respectively; and a is the sJ»oek velocity. A< denotes tlwMuch number. 


3. Sunilarity iittliiiMMW 


The siinilarity tranrfomtations forthe prot^cih Uioler 


hexppnt) 

en»i). 
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4. Numefkiti Sttitttiotii 


For e*lubiun({ the numerical solulionit, it 
the oon-tlimcniiiurutl forms as 


mcoh’ 






The variable /? assumes Ihe pluc.’Jal r^'\te w'esnamttn 

No« ,»ni| Ite ^ Md (l 3 HI?jii«<.;tl» tU,mk^(S):Z ,« ^ 

^ ^‘1 ■■ .: 

1 , 1' '\4 ( 18 ) 

with /I » ^ a remate 


“*■>. ^ _ ' - -: rWl-' 

* ^ ® ^^*^nsioiit^ra<hiatioB parameter. 

.h. Cid.u IIMJ) .,. .,.„di,n»l W ': 

.' . • ' ?, ■ • , V. . 

- t/jfi - F)[(c** '•>Q / Nt*) - FI A,* 

nC(^~- n’ - fp/or '• 


n[(n - Fi’ - (p/c)j 

0 ' « 

idi-y) 

P' « 0[(n - F) F' - F], 

(?' - <Ci,[r{(i, - F)» - lyP/C)} - FW - F>J - 2nu F fV} 

where primes dciioios ditfcrcntiatlon with respa*t to ip 
The appropriate transformed shock conditions arc 

fM 
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umcncul iiUcgitiium i» cumoii uui uniU tlw kiiiemaiw conUi 
•d and il ha. been perfonned on UliS-Sy.iem \m con.puler 
Kanpur by Ihe ^^cll-known RKUS pforanmve for y » 1.4. M 
DO at ■» l. The nature of the lickf variables is illustrated in Tab 
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